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ABSTRACT 

We present the discovery of 11 new double degenerate systems containing extremely low-mass white dwarfs 
(ELM WDs). Our radial velocity observations confirm that all of the targets have orbital periods < 1 day. We 
perform spectroscopic fits and provide a complete set of physical and binary parameters. We review and com¬ 
pare recent evolutionary calculations and estimate that the systematic uncertainty in our mass determinations 
due to differences in the evolutionary models is small (ft; 0.01 Mq). Eive of the new systems will merge due 
to gravitational wave radiation within a Hubble time, bringing the total number of merger systems found in the 
ELM Survey to 38. We examine the ensemble properties of the current sample of ELM WD binaries, including 
the period distribution as a function of effective temperature, and the implications for the future evolution of 
these systems. We also revisit the empirical boundaries of instability strip of ELM WDs and identify new 
pulsating ELM WD candidates. Einally, we consider the kinematic properties of our sample of ELM WDs and 
estimate that a significant fraction of the WDs from the ELM Survey are members of the Galactic halo. 

Keywords: binaries; close - Galaxy: stellar content - gravitational waves - supernovae: general - techniques: 
spectroscopic - white dwarfs 


1. INTRODUCTION 

Short period binary white dwarfs (WDs) are the proposed 
progenitors of transient events such as supernovae la, un- 
derluminous .la, and Ca-rich supernovae, and other exotic 
systems like AM CVn, R Coronae Borea l is (R CrB), and 


single subdwarf B/O stars (IWebbinkI 119841 

Iben & TutukovI 

19841; iBildsten et al.ll2007l; 

Perets et al.ll201C 

;lSolheimll2010l; 

Clavtonll2013HFolevl 12015 

). These binary WDs, including 


the interacting AM CVn systems, are also expected to be ex¬ 
cellent gravitational wave sources and the only known veri¬ 
fication binaries f or the evolved Laser Interfe rometer Space 
Antenna (eLISA, lAmaro-Seoane et aH |2012|) . The initial 
searches for short period doub le WDs have foun d systems 
with periods as short as 1.4 h (iMoran et al.lll997l) . but they 
failed to find a large number of merging binary systems 
(iNaniwotzki et al.ll2007li . 


Low-mass WDs with M ^ 0.4 Mq form » 10% of the 
WD population in the solar neighborhood (iLiebert et alJ 
l2005h, and the majoritv of them are found in binary sys- 


terns ( 

Marsh et alJI 199.51 

Maniwotzki et al.ll2007t IBrown et alJ 

|2mla 

;lDebes et al.ll2015 

). This is expected, as the Galaxy is 


not old enough to produce M < 0.5 Mq WDs through sin¬ 
gle star evolution. The young est WDs in Milky W ay’s globu¬ 
lar clusters have M > 0.5 Mq (iHansen et al.ll2()0% . consistent 
with this explanation. 

Extremely low-mass (ELM) WDs with M < 0.3 Mq pro¬ 
vide a unique opportunity to significantly enlarge the known 
population of merging WDs in the Galaxy. The ELM Survey 
is a targeted search for these WDs in a well defined color and 
tempera ture r ange (Brown et al.l2010l 120121120131: iKilic et aTI 
l2010bl I 2 OI iL 12012 ). This survey has so far identified 56 hi- 


Based on observations obtained at the MMT Observatory, a joint facil¬ 
ity of the Smithsonian Institution and the University of Arizona. 

^ Visiting Astronomer, Kitt Peak National Observatory, National Opti¬ 
cal Astronomy Observatory, which is operated by the Association of Uni¬ 
versities for Research in Astronomy (AURA) under cooperative agreement 
with the National Science Foundation. 


naries, all with P < 1 d, including 33 merger systems and 
four binaries with P < 1 h. The two shortest period systems, 
J0651 and J0935 (WD 0931 +444), will merge in < 1 Myr 
and < 10 Myr, respectively (iBrown et aLll201 Ibl IKilic et~ffl 

120I41 . 

ELM WDs display a v ariety of photometr ic effects, in¬ 
cluding doppler beaming (IShporer et al.l 120101) . tidal distor- 
tion s, pulsations, and eclipses. There are currently six eclips¬ 
ing (Steinfadt et al. I20 id [Brown et all 1201 Ibt IVennes et al 


iParsons et al. 1201 1 : IKilic et al. 2014bt Hallakoun et al 


2011 


2015 ), six pulsating (iHermes et al.ll2013al IKilic et al.ll2015h . 
and eight tidally distorted (IHermes et al.ll2014li ELM WDs 
known. The discovery of new ELM WDs, hence, provides 
new opportunities to improve the mass-radius relation and our 
u nderstanding of the i nterior structure of low-mass WDs. 

IBrown et al.l (120121) describe an efficient way to delineate 
ELM WDs from normal WDs, A stars, and quasars, by target¬ 
ing a color-selected sam ple of g = 15-20 m ag low-mass WD 
candidates (see Eig. 1 of iBrown et al.ll20T2h . Here we extend 
the target selection to the SDSS Data Release 9 area while 
continuing to apply the same color selection criteria and pro¬ 
vide the first discoveries from this extended survey. Three of 
our targets, J11304-3855, J1526H-0543, a ndJ1617-Hl310 were 
previously identified as ELM WDs by IBrown et al.l (120121) 
based on single epoch observations. Here we present follow¬ 
up radial velocity observations of these three targets, as well 
as nine newly identified ELM WDs. 

We have also been obtaining radial-velocity measure¬ 
ments for candidates from other sources including the Large 
Sky Area Multi-Obiect Spectros c opy Telescope (L AMOST, 
iWang et al.ll996tlCui et al.l2012l) . iZhao et ^ (120131) describe 
how a number of WDs were identified from the LAMOST 
data but objects with logg < 7.0 were rejected as WDs. The 
spectra of these discarded objects were kindly provided to 
us by J. K. Zhao (2013, private communication). Based on 
our fits of these spectra, we identified a number of possible 
ELM WD candidates including 103084-5140 and 112494-2626. 
Coincidentally, 112494-2626 lies within the SDSS footprint, 
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Figure 1. ID model fits (red) to the obsei'ved Balmer line profiles (black) for the 11 new ELM WDs binaries. In the right panel, we show fits to 9 ELM WDs, 
with optical spectra obtained at the MMT and FLWO, using pure H model atmospheres and include the lines from H 7 (bottom) to H12 (top). We also show model 
fits for J1054—2121 (left panel, top) whose spectrum was obtained using the Kitt Peak 4 m and includes the lines from H/3 to H9. Finally, fits using a mixed H-He 
model grid are shown for J0756+6704 (left panel, bottom) where we fit the He I A4026 and He I A4471 lines in addition to the Balmer lines. Individual spectral 
lines are offset by a factor of 0.2 for clarity. The best-fit atmospheric parameters are indicated at the bottom of each panel. 


however J0308+5140 does not. 

Section 2 describes our spectroscopic observations. Sec¬ 
tion 3 presents the physical and orbital parameters of the 
eleven new ELM WDs in our survey. The sample charac¬ 
teristics of all 67 binaries in the ELM Survey are discussed in 
Section 4 and we conclude in Section 5. 

2. OBSERVATIONS AND MODELLING 
2.1. Optical Spectroscopy 

We used the 6.5m MM T telescope equipped with the 
Blue Channel spectrograph (iSchmidt et al.l f T98^ . the 200- 
inch Hale telescope equipped with the Double spectrograph 
dOke & Gunnlll982h . the Kitt Peak National Observatory 4m 
telescope equipped wit h the R-C spectrogra ph, and more re¬ 
cently with KOSMOS (iMartini et al.ll20 A . to obtain spec¬ 
troscopy of our eleven targets in several observing runs. We 
operated the Blue Channel, Double, R-C, and KOSMOS spec¬ 
trographs with the 832 line mm“', 1200 line mm“', KPC22B, 
and b2k gratings, providing wavelength coverages of 3650- 
4500 A, 3650-5180 A, 3700-5170 A and 3500-6200 Aand 
spectral resolutions of 1.0 A, 1.7 A, 2.2 A, and 2.0 A, respec¬ 
tively. We also observed two of the targets with g < 17 mag 
in queue scheduled time at t he 1.5 m ELWO tele scope us¬ 
ing the PAST spectrograph (iPabricant et al.l [TOOSh with the 
600 line mm“' grating, providing 3500-5500 A wavelength 
coverage with a spectral resolution of 2.3 A. The Kitt Peak 
and Palomar observations were obtained as part of the NOAO 
programs 2012A-0055, 2012B-0114, 2013A-0276, 2013B- 
0130, 2014A-0189, and 2015A-0082. The MMT and Hale 


observations were obtained at the parallactic angle, but the 
Kitt Peak and ELWO observations were obtained at a fixed 
slit angle. A comparison lamp exposure was obtained with 
every observation . We flux-calibrate u sing blue spectrophoto- 
metric standards (iMassev et al.lll988h . 

We measure r adial velocities using the cross-correlation 
package RVSAO (iKurtz & Minkll998 ). Given the differences 
in wavelength coverage and resolution of the five instruments 
used in our work, we cross-correlate the observed spectra with 
the best-fit model WD templates (see §2.2) at the appropri¬ 
ate resolution for each spectrum. The average precision of 
our measurements ranges from about 10 km s“' to 20 km 
s~*. We compute best-fi t orbital elements using the code of 
iKenvon & Garcial (119861) . which weights each velocity mea¬ 
surement by its associated error. We perform a Monte Carlo 
anal ysis to verify the u ncertainties in the orbital parameters 
(see lBrown et al.l[2012h . 

2.2. WD Model Atmosphere Analysis 

Pigure[T]presents the Balmer line profile fits for our 11 new 
ELM WDs. We use a pure-hydrogen model atmosphere grid 
covering Teff = 4000 - 35,000 K and l ogg = 4.5 - 9.5 and the 
spectro scopic technique described in iGianninas et akl (1201IL 
l2014ai and references therein) to ht the Balmer line profiles 
of each ta rget. The models include the Stark broadening pro¬ 
files from iTremblav & Bergeronl (120091) . We ht all of the vis¬ 
ible Balmer lines up to H12 in the MMT and Hale data, but 
we only ht H/3 through H9 for the Kitt Peak data due to the 
decrease in sensitivity below 3800 A. The results from differ¬ 
ent telescopes agree within the errors, but we adopt the values 
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Table 1 

ELM WD Physical Parameters 


SDSS 

R.A. 

Decl. 

Teff 

logs 

Ml 

^0 

M, 

d 

■^0001 


(h:m:s) 

((d:m:s) 

(K) 

(cm s“^) 

(Mq) 

(mag) 

(mag) 

(kpc) 

(Gyr) 

J0308+5140 

03:08:18.19 

+51:40:11.52 

8380 ± 140 

5.51 ±0.10 

0.151 ±0.024 

13.049 ± 0.010 

8.04 ± 0.53 

0.100 ±0.025 

1.128 ±0.358 

J0756+6704 

07:56:10.71 

+67:04:24.69 

11640 ±250 

4.90 ±0.14 

0.181 ± 0.011 

16.233 ± 0.020 

5.22 ± 0.46 

1.597 ± 0.343 

0.189 ±0.078 

J0837+6648 

08:37:08.51 

+66:48:37.12 

11400 ±240 

6.31 ± 0.05 

0.181 ± 0.011 

17.844 ± 0.019 

8.83 ± 0.25 

0.634 ± 0.072 

1.157 ± 0.123 

J1054-2121 

10:54:35.78 

-21:21:55.94 

9210 ± 140 

6.14 ±0.13 

0.168 ±0.011 

18.487 ± 0.013 

9.11 ±0.49 

0.751 ±0.169 

2.649 ± 0.256 

J1108+1512 

11:08:15.51 

+ 15:12:46.74 

8700 ± 130 

6.23 ± 0.06 

0.167 ± 0.010 

18.830 ± 0.017 

9.61 ±0.28 

0.698 ± 0.090 

3.602 ± 0.505 

J1130+3855 

11:30:17.45 

+38:55:50.11 

11430 ± 190 

7.12 ±0.06 

0.286 ±0.018 

19.446 ±0.021 

10.34 ±0.26 

0.662 ± 0.079 

0.287 ± 0.055 

J1249+2626 

12:49:43.57 

+26:26:04.22 

10120 ± 160 

5.72 ± 0.05 

0.161 ±0.010 

16.566 ±0.015 

7.76 ± 0.24 

0.576 ± 0.064 

1.541 ±0.153 

J1449+1717 

14:49:57.15 

+ 17:17:29.33 

9700 ± 150 

6.08 ± 0.05 

0.168 ± 0.010 

17.620 ± 0.018 

8.76 ± 0.25 

0.591 ± 0.068 

2.152 ±0.128 

J1526+0543 

15:26:51.57 

+05:43:35.31 

10290 ± 190 

5.69 ± 0.05 

0.162 ±0.010 

18.754 ± 0.020 

7.63 ± 0.25 

1.677 ± 0.193 

1.413 ±0.150 

J16I7+1310 

16:17:22.51 

+ 13:10:18.87 

10510 ± 170 

6.07 ± 0.05 

0.172 ±0.010 

18.602 ± 0.015 

8.47 ± 0.24 

1.062 ±0.117 

1.626 ± 0.099 

J2151+1614 

21:51:59.21 

+ 16:14:48.72 

12300 ± 230 

6.24 ± 0.06 

0.176 ± 0.010 

16.454 ± 0.014 

8.51 ±0.24 

0.387 ± 0.044 

0.882 ± 0.098 


Notes. The listed Teff and logg values are corrected for 3D effects according to^embla^^h 



g - r (mag) 


Figure 2. Color-color diagram of the WDs in the ELM Survey (blue circles), 
including 10 of the 11 systems identified in this paper (red triangles). Ma¬ 
genta lines show synthetic colors of WD model sequences with Teff = 30,000 
- 7500 K and logg = 5, 6, 7, and 8. The two shortest period detached binary 
WD systems currently known, J0651 and J0935 (WD 0931 +444), are labeled. 


from the MMT data (when available) since those spectra have 
the highest signal-to-noise ratio and resolution, and include 
all of the higher order Baimer lines. 

J0756+6704 also displays the He I A4026 and He I A4471 
lines in its optical spectrum. In this case,we use a grid of 
mixed H-He model atmospheres to fit T^ff, logg, as well as the 
helium abundance. The measured abundance of log (He/H) 
= -0.60 ± 0.38 coupled with the rather low surface gravity 
of logg= 4.89 m akes it very similar to the mixed H-He ELM 
WDs analyzed in iGianninas et al.l(l2014al) and is potentially a 
sign of a recent shell flash. 

We also note the continuing presence of the Ca II A3934 
line in the blue wing of He for all ELM WDs with logg 
<6.0 (iBrown et al.ll2013t iGianninas et al.ll2014ah including 
J0308H-5140, J0756H-6704, J1249 h-2626, J1526h- 0543 as well 


as J1054-2121. Eurthermore, in both our spectrum and the 
SDSS spectrum of J0756H-6704, the Mg II A5173-5184 dou¬ 
blet is clearly visible. The Na D lines are also definitively 
detected in the SDSS spectrum of J0756 and the LAMOST 
spectrum of J0308-I-5140. These features mark J0756-I-6704 
and J0308-I-5140 as mo st intriguing systems worthy of a more 
detail ed analysis (e.g. IGianninas et akl l2014at iHermes et al.l 
l2014ft . In all cases, the spectral range where the Ca line is 
present is excluded from both the normalization and fitting 
procedures. 

Table [T] summarizes the measured atmospheric parameters, 
Teff and logg. Note that the values listed in Table [T] have 
been corrected using the new 3D correction functions from 
iTremblav et al.l (1201 5h and thus differ from those displayed 
in Eigure [T] The corrections have been applied to all targets 
with Teff ^ 12,000 K which includes all 11 new ELM WDs 
save J2151H-1614. 

Based on the r e cent e volutionary sequences for ELM WDs 
by lAlthaus et al.l (120131) . we estimate masses (Mi), absolute 
magnitudes (Mg), and WD cooling ages (Tcooi) for each object. 
Given Mg and the extinction corrected SDSS g-band magni¬ 
tude (go). However, since J0308-I-5140 is located outside the 
SDSS footprint, for this object we estimate the g-band magni¬ 
tude from The Eou rth US Naval Observat ory CCD Astrograph 
Catalog (UCAC4, IZacharias et al.ll2M^ . We also derive dis¬ 
tances (d) to each object. These are also included in Table 
□ All 11 targets have Teff sa 8000 - 12,000 K, logg = 4.9 - 
7.1, M sa 0.15 - 0.29 Mq and they are located at d sa 0.10 - 
1.70 kpc from the Sun. __ 

Contrary to previous analyses (IGianninas et al.ll2014allbh . 
we do not adopt a systematic uncertainty of 0.02 Mq for our 
determination of Mi. Instead, we propagate th e uncertainties 
in Tpff and logg through our interpolation of the lAlthaus et al.l 
(120131) model grid after which we add, in qu adratu re, the sys¬ 
tematic uncertainty of 0.01 Mq (see Section l4.2.11 i. 

Eigure |2] shows a color-color diagram of all of the ELM 
WDs identified in the ELM Survey, including 10 of the 11 sys¬ 
tems described in this paper. The latter systems were selected 
to have g-r = -0.2 t o 0.0 mag based on the color-selection 
of iBrown et al.l (120121) . Hence, it is not surprising that they all 
have temperatures near 10,000 K. This figure demonstrates 
that ELM WDs have colors consistent with the WD model se¬ 
quences, except J0935 (WD 0 931-1-444), whose p hotometry is 
contaminated by an M dwarf (iKilic et al.ll2014al) . 
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Figure 3. Periodograms for the 11 new ELM WD binaries. The best orbital 
periods have the smallest values; some binaries are well constrained and 
some have period aliases. 


2.2.1. J0308+5140 

Fitting the optical spectrum of J0308+5140 required an 
additional treatment. It is clear from the spectra that 
J0308+5140 suffers from a considerable amount of redden¬ 
ing, consistent with its Galactic latitude of b = -5.67°. Ac¬ 
cording to the NASA/IPAC Infrared Science Archive (IRSA) 
dust mapfl the reddening in the direction of J0308-I-5140 is 
E{B-V) = 0.795. Adopt ing this value f or the reddening, we 
follow the prescription of lSeatorJ (Il979h to de-redden the ob¬ 
served spectrum of J0308-I-5140. The slope of the corrected 
spectrum is in excellent agreement with the slope of the spec¬ 
troscopic solution. Consequently, we adopt an extinction of 
Ag = 3.014 from IRSA to correct the g-band magnitude of 
J0308H-5140fromUCAC4(i.e. g= 16.063)togetgo= 13.049, 
as listed in Table [T] 

3. ORBITAL PARAMETERS 

Figure [3] shows the periodograms for our 11 targets. All 
of our targets have a well defined minimum in the peri¬ 
odograms and the period is constrained to P < 1 d in all cases. 
Furthermore, to obtain a p-value of p = 0.01 with four degrees 
of freedom requires a Ax^ > 13.3, with respect to the mini¬ 
mum x^ (iPress et al.lll986^ . This criterion is met by all the 
new ELM WD binaries. Thus, the orbital periods we deter¬ 
mine are statistically significant at the 99% confidence level. 

Table |2] and Figure |?] present the best-fit orbital parame¬ 
ters and the radial velocity curves, respectively. Table |2] lists 
the orbital period (P), radial velocity semi-amplitude (K), sys¬ 
temic velocity ( 7 ), the mass function, the minimum secondary 
mass (M 2 ) assuming / = 90°, the secondary mass assuming 
i = 60°, the maximum gravitational wave merger time (Tmerge), 
the orbital separation (a) and finally, the gravitational wave 
strain (log/i). 

Our targets show radial velocity variations with orbital pe- 

* http://irsa.ipac.caltech.edu/applications/DUST/ 
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Figure 4. Obseiwed velocities phased to the best-fit orbits (red) for the 11 
new ELM WD binaries. 


riods of P Ki 3-19 h and semi-amplitudes of /T 80 to Ri 
285 km s“'. The systemic velocities range from 7 w -100 to 
Ri 100 km s“'. We note that the systemic velocities do not 
include the corrections for the WDs’ gravitational redshift, 
which is a few km s“'. 

The minimum mass companions to these 11 ELM WDs 
range from 0.16 Mq to 0.82 Mq. Based on the mass func¬ 
tion, the probability of a 1.4-3.0 Mq neutron star companion 
ranges from 2 % to a maximum of 20 % for the companion of 
J0756-I-6704. These percentages are quite low and it is un¬ 
likely that any of these 11 new ELM WDs have neutron star 
comp anions. These probab ilitie s agree with th e recent stud¬ 
ies o f lAndrews et al.l (l2014ll and iBoffinI (12015h whose results 
imply that neutron star companions to ELM WDs should be 
rare. 

It is also unlikely that the total masses for any of these new 
systems is above the Chandrasekhar mass limit. Assuming 
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Table 2 

ELM WD Binary Parameters 


SDSS 

P 

(days) 

K 

(km ) 

7 

(kms-*) 

Mass Function 
(Mq) 

M2 

(Mq) 

3^2,1=60° 

(Mq) 

Emerge 

(Gyr) 

a 

(Re) 

\ogh 

J0308+5I40 

0.80590 ± 0.00038 

78.9 ± 2.7 

-56.8 ± 1.8 

0.041 ± 0.004 

S; 0.16 ±0.02 

0.20 ± 0.03 


2.46 ±0.12 

-22.57 

J0756+6704 

0.61781 ±0.00002 

204.2 ± 1.6 

-98.7 ± 1.3 

0.545 ±0.013 

S; 0.82 ± 0.03 

1.13 ±0.04 


3.05 ± 0.04 

-23.05 

J0837+6648 

0.46329 ± 0.00005 

150.3 ± 3.0 

-22.4 ± 2.3 

0.163 ±0.010 

^ 0.37 ± 0.02 

0.48 ± 0.03 


2.06 ± 0.04 

-22.85 

J1054-2121 

0.10439 ±0.00655 

261.1 ± 7.1 

94.3 ± 4.5 

0.193 ±0.028 

Ss 0.39 ± 0.05 

0.52 ± 0.06 

sS 1.452 

0.77 ± 0.06 

-22.49 

J1108+1512 

0.12310 ±0.00867 

256.2 ± 3.7 

44.5 ± 2.4 

0.214 ± 0.024 

S; 0.42 ± 0.04 

0.56 ± 0.05 

sj 2.143 

0.87 ± 0.07 

-22.48 

J1130+3855 

0.15652 ±0.00001 

284.0 ± 4.9 

24.5 ± 3.6 

0.371 ±0.019 

S; 0.72 ± 0.04 

0.96 ± 0.05 

s; 1.652 

1.23 ±0.02 

-22.14 

J1249+2626 

0.22906 ±0.00112 

191.6 ±3.9 

-15.5 ±2.4 

0.167 ±0.011 

> 0.35 ± 0.02 

0.47 ± 0.03 

s; 13.242 

1.26 ±0.03 

-22.65 

J1449+1717 

0.29075 ± 0.00001 

228.5 ± 3.2 

3.9 ±3.0 

0.359 ±0.015 

> 0.59 ± 0.03 

0.81 ±0.04 


1.69 ±0.03 

-22.54 

J1526+0543 

0.25039 ± 0.00001 

231.9 ±2.3 

-39.8 ± 1.8 

0.324 ± 0.010 

^ 0.54 ± 0.02 

0.74 ± 0.02 

s; 12.087 

1.49 ±0.02 

-22.99 

J1617+1310 

0.41124 ±0.00086 

210.1 ±2.8 

7.7 ± 2.0 

0.395 ± 0.017 

S; 0.64 ± 0.03 

0.87 ± 0.04 


2.17 ±0.04 

-22.86 

J2151+1614 

0.59152 ± 0.00008 

163.3 ±3.1 

-67.0 ± 2.3 

0.267 ± 0.015 

> 0.49 ± 0.03 

0.66 ± 0.04 


2.59 ± 0.05 

-22.61 



Figure 5. Location in the Teff - log^ plane of the 11 new ELM WDs pre¬ 
sented in this paper (red triangles) along with the complete ELM Survey 

S ’; (blue circles). Als o plotted are evoluti onary tracks from lAlthaus et all 
(solid lines) and llstrate et alj 420141) (dashed lines) for 0.165 Mq 
(cyan), 0.187 Mq (magenta), and 0.239 Mq (black). 

an average inclination of i = 60°, the binaries with the high¬ 
est total system mass are J0756H-6704 and J1130H-3855 with 
Mtot =1.31 Mq and 1.25 Mq, respectively. 

Our sample also includes five new merger systems 
(J1054-2121, J1108H-1512, J1130 h-3855, J1249h- 2626, 
J1526H-0543). The quickest to merge will be J1054-2121 
with Tmerge = 1.452 Gyr. However, the gravitational wave 
strains of these systems ar e not strong enough to b e detected 
by eLISA (see Figure 12 of iGianninas et al.ll2014ah . 

4. DISCUSSION 

4.1. Eleven New Binary WD Systems 

With these 11 new discoveries, and including the p ul- 
sating ELM WDs published by [Hermes et al.l (l2013bllah . 
and the discoveries presented in lKilicetal.1 (l2014albl) and 


^Althaus (^q) 

Figu re 6. Difference in t he primary m ass as determined u sing the models 
from lAlthaus et all 420131 MAithaus) and llstrate et alj 420141 Mistrate) plotted 
as a function of Teff (top) and MAithaus (bottom). The dotted line denotes the 
average difference. 


iGianninas et al.l (I2014bl) . the ELM Survey has found a total of 
73 ELM WDs of which 67 are in detached, double-degenerate 
binaries; 38 of the binaries will merge within a Hubble time. 
Counting only six systems where no significant radial veloc¬ 
ity variability has been detected, 92% of the ELM Survey 
targets are formally members of compact (P < 1 d) binary 
systems. In addition to the results listed in Tables[T]and|3 we 
provide, for completeness, the physical and binary parameters 
for the entire ELM Survey sample in Tables |4] and |5] of the 
Appendix. Note that for ELM WDs with Teff < 12,000 K, 
the atmospheric parameters have been adjusted using the 3D 
corrections from iTremblav et Tl (120151) . 
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Figure 7. Cumulative distribution as a function of T^ff for ELM WDs with 
P < 0.1 d (solid line) and P > 0.1 d (dotted line). 

4.2. The ELM WD Sample 
4.2.1. Comparison of Evolutionary Models 

In Table [T] we adopt masses derive d from the evolution - 
ary sequences for He-core WDs of lAlthaus ^ ID (Eoil. 
However, new evolutionary sequences appropr iate for He- 
core E LM WDs have recently been published bv llstrate et alJ 
(l2014l) . Due to the importance of obtaining accurate masses 
in order to correctly predict the final merger products of ELM 
WD binaries, it is important to explore these new evolutionary 
calculations. 

One of the major diffiren ces between the lAlthaus et alJ 
dMl and llstrate et alJ (12014|) models are the initial assump¬ 
tions with regards to the proge nitors and companion s of the 
eventual ELM WDs. All of the lAlthaus et al.l d2013l) models 
assume a 1.0 Mq progenitor for the ELM WD and are evolved 
in a binary with a lAM g^ neutron star companion. In contrast, 
the llstrate et’aP (12014ft models assume a range of progenitor 
masses (from 1.2 to 1.6 Mq) as well as a range of neutron 
star masses (from 1.3 to 1.75 Mq) to produce the full range 
of ELM WD masses in their model grid. As we will demon¬ 
strate, both approaches yield similar ELM WD masses. How¬ 
ever, as we mentioned in Section 3, neutron star companions 
are unlikely in most of these systems. Consequently, it would 
be useful to have evolutionary tracks computed for ELM WD 
binaries with massive WD companions instead to provide ad¬ 
ditional insight into the final ELM WD masses. 

In Eigure|5] we plot the location of all ELM WDs from the 
ELM Survey sample in the Teff - logg plane. I n addition, we 
plot a selection of evolut ionary tracks from both lAlthaus et al.l 
loll and llstrate et al.l (120141) . In order to make the most di¬ 
rect comparison possible, we have chosen a set of three repre¬ 
sentative evolutionary tracks corresponding to final ELM WD 
masses of 0.165 Mq, 0.187 Mq, and 0.239 Mq. 

The results in Eigure|5]demonstrate that there are significant 
differences between the two set s of evolutionary tracks. Most 
notably, the llstrate et'd] (12014ft models predict considerably 
fewer shell flashes than the lAlthaus et al.l (120131) tracks. This 
is most clearly seen in the 0.187 Mq tracks (magenta) where 
the Istrate model predicts no shell flashes while the equivalent 


Figure 8. Plot of the minimum total mass as a function of the maximum 
merger time for the entire ELM Survey sample (blue circles) with the 11 new 
binaries from this paper (red triangles). 

Althaus model undergoes a series of shell flashes. Indeed, 
the Althaus models assume Z =0.01 while the Istrate models 
assume Z =0.02. Since the maximum m ass for which flashe s 
occur increases with lower metallicity (iNelson et al.l 120041) . 
it is not surprising that the Althaus models predict a greater 
number of shell flashes and for a broader range of masses. 

Despite these significant differences, the most important 
question is whether or not the models yield the same masses 
for a given Teff and l ogg. We have chos en a set of nine mod¬ 
els from the grid of llstrate et alJ (12014ft to compute masses. 
To greatly simplify the calculation, we consider only the final 
cooling portion of the evolutionary tracks. In Eigure|6]we plot 
the difference in masses derived from the two model grids as 
function of and MAithaus- Eigure |6] demonstrates that the 
Althaus & Istrate models yield very similar masses. The aver¬ 
age difference is Ri 0.01 Mq with an overall scatter of roughly 
the same magnitude. Consequently, we adopt an uncertainty 
of 0.01 Mq as the systematic error due to uncertainties in the 
evolutionary models. This is equivalent to an error of Ri 5% in 
mass. In contrast, the mean error in mass due to uncertainties 
in the atmospheric parameters, Teff and logg, is Ri 0.006 Mq 
which is « 3%. Therefore, the error in our mass estimates is 
dominated by the systematic uncertainties from the evolution¬ 
ary models. 

4.2.2. The Period Distribution of Binary WDs 

In Eigure |2] we plot the cumulative distributions, as a func¬ 
tion of Teff, for the 18 systems with P < 0.1 d and the 49 
systems with P > 0.1 d. We see that the majority of systems 
with P > 0.1 d are found with Teff < 15,000 K. 

To determine if the two distributions are statistically inde¬ 
pendent, we perform the two-sample Kolmogorov-Smirnov 
(K-S) test. The K-S statistic, is defined as the largest 
difference between the cumulative distributions functions and 
can be expressed as 

Ai.«' =SUp\Fy„(x)-F2,n'{x)\ 

X 

where Ei_„(.r) andT 2 „'(x) represent the P > 0.1 d andP < 0.1 d 
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Figure 9. Plot of M\ vs. M 2 for the entire ELM sample. Black points repre¬ 
sent Mo.mm while the red triangles assume M 2 ,,= 60 ° • The majority of the ELM 
WD binaries lie in the region between the regions of stable and unstable mass 
transfer. 


distribtuions, respectively. Eor the given distributions in Eig- 
ure I 2 I we determine a K-S statistic of „/ = 0.5590 and a 
p-value of p = 0.0002. Based on this result, we can reject 
the null hypothesis at the a = 0.01 significance level. How¬ 
ever, the K-S test is not always sensitive enough to determine 
if two distributions are indeed independent. This is because, 
by definition, the distributions converge to 0 and 1 at the ends. 
The Anderson-Darling (A-D) test is more useful in such situa¬ 
tions. We perform the A-D test and obtain a standardized test 
statistic of r = 6.9807 and p = 0.0006. Based on these values, 
we once again reject the null hypothesis at the a = 0.01 sig¬ 
nificance level. The results of these statistical tests strongly 
suggests that the two samples are not drawn from the same 
parent distribution. 

The fact that the two samples appear to originate from sep¬ 
arate parent distributions can be understood intuitively. The 
shorter period systems (E < 0.1 d) will merge before they 
have a chance to cool to Teff < 15,000 K. On the other hand, 
the longer period binaries (P > 0.1 d) have maximum merger 
times which are considerably longer and therefore they can 
survive long enough to cool below Teff ^ 15,000 K. 

4.2.3. The Future: Merger Products 

In Eigure 0 we plot the total minimum system mass (as¬ 
suming i = 90°) as a function of the maximum merger time of 
the system. The maximum gravitational wave merger time is 
given by 


Emerge,max 


M.\M.2 


P&/3 X lO-^Gyr 


( 1 ) 


where Mi is the mass of the ELM WD in Mq, M 2 is the min- 
imum companion mass i n Mq, and the period P is in hours 
(iLandau & Lifshit3ll958h . 

Our new sample of 11 ELM WD binaries has yielded 
five new systems (J1054-2121, J1108H-1512, J1130H-3855, 
J1249H-2626, J1526H-0543) which are expected to merge 
within a Hubble time. Thus, the ELM Survey has now identi¬ 
fied a total of 38 merger systems. 

10651 and WD 0931H-444 remain the two shortest period 
systems, expected to merge in « 1 Myr and w 10 Myr, re¬ 
spectively. The most massive merger systems J1741H-6526, 



Figure 10. Region of the Teff - logg plane containing the canonical ZZ Ceti 
instability strip (lower left) and the six currently known ELM pulsators (upper 
right). All pulsators are identified as open diamonds, whereas WDs that have 
been confirmed as photometrically constant are repre sented as filled dots. Th e 
solid lines denote the blue and red boundaiies from iGianninas et al.ll2Q14af) . 
The blue and red dashed lines indicate the tentative new boundaries which 
take into account the 3D corrected values of Tgff and logg. Error bars denote 
WDs from the ELM Survey sample which have not yet been investigated 
for photometric variability. Red triangles indicate the 11 new ELM WDs 
binaries. 


10751-0141, near the top of Eigure[8] represent potential Type 
la progenitors. However, the total mass of the system is not 
the o nly parameter that d etermines the ultimate fate of the bi¬ 
nary. iMarsh et'aH (120041) showed that systems with high mass 
ratios will in fact lead to stable mass tran sfer binaries. In Eig- 
ure |9] we plot the stability diagram from iMarsh et~an (l2004l) 
and also plot the values of Mi and M 2 for the entire ELM Sur¬ 
vey sample for both the minimum and most likely (i.e. for 
i = 60°) companion mass. Note that in all cases, we choose 
M2 > Mj^ _ 

As iKilic et aTI (12014bh have shown, 11741-(-6526 and 
10751-0141 are, given the uncertainties in M 2 , located in the 
region of the stability diagram which predicts that they will 
eventually become stable mass transfer systems. Thus, they 
represent the progenitors of future AM CVn systems. How¬ 
ever, Eigure |9] also shows that two more ELM WD binaries lie 
this same region, namely 10811-(-0225 and 12132-(-0754. This 
brings the total to four unambiguous identifications of AM 
CVn progenitor systems. In addition, a number of systems 
do lie in the the unstable regime, which predicts an eventual 
merger. However, the total mass for all of these merger sys¬ 
tems is significantly below the Chandrasekhar mass. 

Eigure |9] also clearly shows that most ELM WD binaries 
are found in the regime between the stable and unstable re¬ 
gions. Indeed, there is still much uncertainty and debate 
over the eventual fate of compact d ouble WD binaries. The 
recent calculations of IShenI (120 15h predict that most dou- 
ble WD bina r ies wi ll eventually merge. On the other hand, 
iKremer et alJ (1201 5h predict that a majority of systems will 
undergo stable mass transfer. Once the ELM Survey is com- 
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plet ed, a more rigorous examination of both the IShenI (12015h 
and iKremer et alJ (1201 5h models will be possible by compar¬ 
ing these merging ELM WD binaries with the known popula¬ 
tions of various merger products including AM CVn systems, 
R CrB stars and Type la supernovae. 

4.2.4. Instability Strip 

The disco very of the fir s t five p ulsating ELM WDs were 
reported in [Hermes etalJ (l2013ah and references therein. 
Recently, iKilic et al.l (120151) discovered the sixth pulsating 
ELM WD which also happens to be th e companion to 
PSR J1738-(-0333 (lAntoniadis et al.l l2012^ . These pulsat¬ 
ing ELM WDs all have = 8000 - 10,000 K and logg 
= 6-7 (see the update d instability strip in Eigure 13 
of iGianninas et alJ 12014ah . Eurthermore, the location of 
this newly defined instability strip is in good agreement 
with the predictions of pulsation models (IVan Grootel et alJ 
l2013h . In Eigure [TO] we plot an updated version of the 
figure from Gianninas et al. including the 11 new ELM 
WDs from this paper and t he canonical ZZ Ceti instability 
strip (IGianninas et al.ll201 Ih includin g several new discover¬ 
ies presented by iGreen et al.l (l2015l) . One important differ¬ 
ence here is that the adopted and logg values for all the 
WDs in the figure have be en corrected fo r 3D e ffects us¬ 
ing the functions defined in [Tremblay et al.| (| 2013h for WDs 
with logg > 7.0 and iTremblav et al.l (1201 5h for the ELM 
WD regime. In particular, the corrected atmospheric pa¬ 
rameters for PSR J1738H-0333 are = 8910 ± 150 K and 
logg = 6.30 ±0.10. 

Our best-fit temperatures and surface gravities place five 
of our new ELM WDs within (J1054-2121, J1108H-1512) or 
near the blue edge (J0308H-5140, J1130H-3855, J1449H-1717) 
of the instability strip as delineated by our new boundaries 
defined by Equations |2] and |2 respectively. Eollow-up high¬ 
speed photometric observations of these five targets will be 
useful to search for pulsations in these ELM WDs and are un¬ 
der way with the McDonald 2.1m telescope (K. Bell 2015, 
private communication). 

(logg)biue = 5.96923 X 10-4(7;ff)biue + 0.52431 (2) 

(logg)red = 8.06630 X 10-4(reff),ed-0.53039 (3) 

4.2.5. Disk vs. Halo Membership 

In Table[3]we list the (U,V, W) space velocities for 49 of the 
73 WDs from the ELM Survey for which proper motions were 
available from the SDSSh-USNO-B catalog of iMunn et ^ 
(l2004ft . In addition, the proper motion for J0345-(-1748 
(NETT 11748) is taken fromiKawka & Venn^ (|2009|) and its 
systemic velocity from iKilic et al.l (l2010al) . Einallv we adopt 
the proper motion from UCAC4 for J0308H-5140. 

We computed the (t/, V, W) space velocities and their 
associated uncertainti e s acc ording to the prescription of 
iJohnson & SoderblomI (Il987h combining the proper motions 
with our determinations of the systemic velocity (see Table|2]l. 
It is important to note that since we do not have parallax mea¬ 
surements for our targets, we adopt our spectroscopically de¬ 
termined distances as a proxy for the parallax. In addition, 
we correct for the motion of the local standard of rest using 
(U, V, IV)(T| = ( 11.10, 12.24, 7.25) km s“' as determined by 
ISchonrich et al.l (120101) . 

In Eigure fin we plot the W vs. V and U vs. V velocity dis¬ 
tributions. In addition, we plot ellipsoids representing Galac- 


Table 3 

Space Velocities of ELM WDs 


SDSS 

U 

(Ion s-') 

y 

(km s“ 

') 

ly 

(Icms-') 

f^m,disk 

^m.halo 

J0022-1014 

82 ±23 

-40 ± 

21 

27 ± 8 

1.944 

1.505 

J0056-06I1 

-4± 7 

-42 ± 

9 

-11 ± 3 

0.520 

1.376 

J0106-1000 

-126 ±51 

-236 ± 

55 

-16± 16 

5.173 

1.122 

JOl12+1835 

67 ± 8 

-104 ± 

9 

52 ± 8 

2.677 

1.055 

JO152+0749 

82 ± 12 

-117 ± 

19 

-15 ± 13 

2.597 

0.812 

J0345+1748 

-147 ± 6 

-197 ± 

25 

-28 ± 4 

4.889 

1.194 

J0651+2844 

- 8 ± 4 

2 ± 

15 

-9± 15 

0.545 

1.795 

J0745+1949 

-67 ± 4 

-10 ± 

4 

23 ± 4 

1.729 

1.797 

J0755+4800 

-23 ± 3 

25 ± 

2 

31 ± 3 

1.440 

2.053 

J0815+2309 

153 ± 19 

-275 ± 

38 

-6 ±25 

6.037 

1.270 

J0822+2753 

64 ± 6 

-20 ± 

9 

-24 ± 7 

1.434 

1.622 

J0825+1152 

-71 ± 18 

54 ± 

24 

-13 ± 24 

2.217 

2.359 

J0849+0445 

-20 ± 15 

-20 ± 

16 

25 ± 17 

0.957 

1.630 

J0900+0234 

- 8 ± 14 

-52 ± 

13 

38 ± 14 

1.363 

1.368 

J0917+4638 

-25 ± 22 

19 ± 

32 

9 ±23 

1.067 

1.973 

J0923+3028 

11 ± 3 

-16 ± 

4 

2± 3 

0.231 

1.618 

J1046-0153 

10 ± 6 

-4± 

7 

-51 ± 6 

1.427 

1.797 

J1053+5200 

-112 ± 22 

-201 ± 

29 

15 ± 13 

4.440 

0.947 

J1056+6536 

20 ±22 

-27 ± 

24 

20 ± 18 

0.769 

1.532 

J1104+0918 

-23 ± 3 

-84 ± 

6 

31 ± 5 

1.716 

1.080 

J1141+3850 

280 ± 91 

-345 ± 102 

106 ± 32 

9.378 

2.664 

J1234-0228 

-2 ± 12 

-70 ± 

12 

25 ± 6 

1.292 

1.157 

J1238+1946 

64 ±47 

-652 ± 

65 

-92 ± 9 

12.952 

4.492 

J1422+4352 

-48 ± 39 

-72 ± 

35 

-165 ± 18 

4.879 

2.071 

J1436+5010 

43 ± 13 

9± 

12 

-26 ± 8 

1.228 

1.871 

J1443+1509 

-33 ± 9 

-167 ± 

22 

-141 ± 5 

4.994 

1.505 

J1448+1342 

6 ± 10 

2 ± 

11 

-21 ± 7 

0.687 

1.795 

J1512+2615 

-5 ± 11 

-60 ± 

13 

-5± 7 

0.824 

1.209 

J1518+0658 

-38 ± 4 

-21 ± 

5 

-35 ± 3 

1.288 

1.645 

J1538+0252 

-179 ± 17 

-126 ± 

24 

-1 ± 19 

4.496 

1.503 

J1557+2823 

26 ± 3 

5 ± 

3 

17 ± 3 

0.889 

1.824 

J1625+3632 

56 ± 57 

-136 ± 

49 

-22 ± 44 

2.639 

0.581 

J1630+2712 

36 ±36 

-91 ± 

34 

-106 ± 30 

3.332 

1.411 

J1630+4233 

34 ± 13 

-8± 

10 

-10 ± 10 

0.727 

1.691 

J1741+6526 

39 ± 14 

-55 ± 

8 

-14 ± 12 

1.072 

1.261 

J1840+6423 

118 ± 18 

-61 ± 

6 

2 ± 12 

2.617 

1.389 

J2132+0754 

48 ± 14 

-32 ± 

10 

- 6 ± 13 

0.993 

1.474 

J2228+3623 

20 ± 5 

-37 ± 

5 

26 ± 5 

0.947 

1.457 

J2236+2232 

-328 ± 33 

-181 ± 

11 

-9± 13 

7.903 

2.447 

J2338-2052 

18 ± 22 

11 ± 

22 

3± 8 

0.702 

1.868 

J2345-0102 

343 ± 35 

-206 ± 

26 

69 ± 15 

8.507 

2.450 

This paper 

J0308+5140 

55 ± 1 

-24 ± 

2 

12 ± 2 

1.189 

1.574 

J0756+6704 

17 ± 21 

-227 ± 

48 

29 ±24 

4.244 

0.523 

J0837+6648 

4± 7 

-53 ± 

10 

8 ± 7 

0.727 

1.276 

J1054-2121 

8 ± 11 

-64 ± 

7 

63 ± 9 

2.080 

1.369 

J1108+1512 

-29± 11 

-40 ± 

11 

24 ± 6 

1.121 

1.454 

J1130+3855 

-85 ± 15 

-29 ± 

12 

-1 ± 6 

1.951 

1.658 

J1526+0543 

21 ±21 

2 ± 

27 

-54 ± 19 

1.575 

1.855 

J2151+1614 

-53 ± 7 

-A4 ± 

3 

14 ± 5 

1.405 

1.452 


tic thick-disk and stellar halo populations. We adopt aver¬ 
age velocities of i{U), (V), (W)) = (4, -20, -3) km s“' and 
(17, -187, -5) km s“', and velocity dispersions of (au, cry, 
aw) = (46, 50, 35) km s”' and (1 41, 106, 94) km s"^ for the 
thick-disk and halo, respectively (IChiba & Beersll2000l) . 

Treating the ellipsoids for the thick-disk and halo as mul- 


























































THE ELM SURVEY VI 


9 




V (km s“0 

Figure 11. W vs. U (top) and U vs. V (bottom) velocity distributions for 49 
ELM WD binaries including nine from this paper. Blue circles denote ELM 
WDs with a Mahalanobis distance suggesting halo membership while the red 
triangles represent ELM WDs with kinematics consistent with the Galactic 
thick-disk. The ellipsoids denote the i-cr (dotted) and 2-cr (dashed) contours 
for Galactic thick-disk and stellar halo populations. 


tivariate normal distributions, we compute the Mahalanobis 
distance, as defined in Equation (|4]l, for each ELM WD. This 
measures the distance from the center of the distributions in 
units of standard deviations and provides a quantitative indi¬ 
cator of group membership. 




i{u-{u)f (y-(y))2 (w-{w))^ 


(4) 


The Mahalanobis distances with respect to the thick-disk and 
halo distributions. Dm,disk and Dm,halo, are listed in Table [3 
Based on these values, 27 ELM WDs have distances consis¬ 
tent with disk membership while the remaining 22 WDs have 
distances consistent with halo membership. J0756-I-6704, the 


ELM WD with the most negative systemic velocity among 
the 11 new ELM WD binaries, numbers among the likely 
halo members. Unsurprisingly, LP 400-22 (J2236-I-2232) also 
stands out in the bottom panel of EigurefTTIwith U = -332 km 
s~*. The case of LP 400-22 is well documen ted (iKawka et al.l 
120061 IVennes et al.l 120091 iKilic et aTI 1201 3h and the results 
here serve only to further confirm that LP 400-22 is indeed 
a unique member of the halo. We also find that NETT 11748 
is a halo objec t , in agr eement with iKawka & VennesI (l2009h 
and iKilic et akl (1201 Oah . 

The breakdown of halo versus thick-disk members suggests 
that Ri 40% of the WDs from the ELM Survey are members 
of the Galactic halo. However, this result is based solely on 
a kinematic metric. Additional considerations, such as the 
location of the ELM WDs versus the scale height of the disk, 
may provide additional information not considered here. As 
such, we suggest that our result represents an upper limit for 
the percentage of halo objects in the ELM Survey. 

Having said that, the fact that our sample includes a sig¬ 
nificant fraction of halo stars should not be unexpected. This 
result is most likely a logical consequence of our candidates 
being largely selected from the SDSS which is a high Galac¬ 
tic latitude survey. Given the luminosity of our targets, our 
apparent magnitude range samples white dwarfs in the range 
04 < |Z| < 2 kpc, where Z denotes the vertical distance from 
the Galactic plane. Given the scale height of the disk and the 
relative normalization of the disk and halo, it is not surprising 
that the ELM Survey contains a significant number of halo 
stars. 


5. CONCLUSIONS 

We present radial velocity measurements and stellar atmo¬ 
sphere fits for 11 new ELM WD binaries with P < 1 d. This 
brings the total of ELM WDs identified by the ELM Survey 
up to 73 demonstrating the continued effectiveness of using 
SDSS colors to identify ELM WD candidates. Based our op¬ 
tical spectra, we perform spectroscopic fits, compute orbital 
solutions, and provide a complete set of physical and binary 
parameters for each system. 

We review the recent evolutionary calculations of 
llstrate et 'm (1201 4h and co mpare them to the existing mod- 
els of lAlthaus et all (l2013h . There are significant differences 
between the two sets of evolutionary models. These can be 
traced back to the different initial assumptions made about 
the nature and composition of the progenitor binaries. We 
can only conclude that there still remains much uncertainty 
regarding the exact evolutionary history of ELM WDs and 
additional modeling is required to form a coherent picture. 
Thankfully, the systematic uncertainty when deriving masses 
from the independent model sets is only ss 0.01 Mq. 

When considering the distribution of orbital periods as a 
function of Teff, we have shown that two distinct populations 
are emerging. The shorter period systems are generally found 
with Teff > 15,000 K while the longer period systems can cool 
to much lower Teff- 

Of the 11 new ELM WD binaries, five will merge within a 
Hubble time. Unfortunately, the final merger product for the 
majority of ELM WD binaries remains uncertain. Nonethe¬ 
less, we have identified four systems which will undergo sta¬ 
ble mass transfer and ultimately become AM CVn systems. 

Based on our updated boundaries of the instability strip, we 
have also identified a number of potentially interesting targets 
for follow-up high-speed photometric observations. These 
ELM WDs have atmospheric parameters within or near the 
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instability strip of pulsating ELM WDs, an extension of the 
canonical ZZ Ceti instability strip. 

Einally, a study of the kinematics of ELM WDs reveals that 
the majority of these binaries are members of the Galactic 
disk. However, a non-negligible fraction of ELM WDs are 
most certainly members of the Galactic halo, a result of our 
reliance on SDSS colors to identify candidates. 
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telescope (RCSpec, KOSMOS), Palomar Hale 5m telescope 
(Double Spec), EWLO 1.5m telescope (EAST) 


APPENDIX 

ELM WD PARAMETERS 

We provide in Tables |4] and |5] a complete listing of the 
physical and binary parameters for the entire ELM Survey 
sample. This includes the 62 ELM WDs from previous pub¬ 
lications as well as the 11 new ELM WDs presented here. 
We note that for ELM WDs with Teff < 12,000 K, the atmo¬ 
spheric pa rameters have been ad justed using the 3D correc¬ 
tions from iTremblav et ^ (l2015h . We also note that we have 
updated our value of the apparent magnitud e for J0345H-1748 
(NETT 11748). In iGianninas et al.l (l2014bl) we had ad opted 
the value of V = 16.5 from iKawka & Vennei (l2009l) . We 
now adopt the value of g = 16.797 from UCAC4. We then 
take A„/E(B — V) = 3.7 9 3 from IRSA and E(B — V) = 0.10 
from iKawka & VennesI (l2009l) to obtain an extinction of 
Ag = 0.3793. Thus, our final extinction-correctedg-band mag¬ 
nitude for NETT 11748 is go = 16.418. 

DATA TABLE 

Table |6] presents our radial velocity measurements for the 
11 new ELM WD binaries presented here. The table columns 
include object name, heliocentric Julian date, heliocentric ra¬ 
dial velocity, and velocity error. 
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Table 4 

ELM WD Physical Parameters 


SDSS 

Tiff 

(K) 

logg 
(cm s“^) 

Ml 

(M©) 

Radius 

(Re) 

go 

(mag) 

(mag) 

d 

(kpc) 

fcool 

(Gyr) 

J0022+0031 

20460 ±310 

7.58 ± 0.05 

0.457 ± 0.016 

0.0182 ± 0.0013 

19.284 ± 0.033 

9.84 ±0.18 

0.774 ± 0.065 

0.215 ±0.128 

J0022-1014 

20730 ± 340 

7.28 ± 0.05 

0.375 ± 0.016 

0.0233 ± 0.0018 

19.581 ±0.031 

9.28 ± 0.20 

1.151 ± 0.107 

0.042 ± 0.021 

J0056-0611 

12240 ± 180 

6.18 ±0.04 

0.174 ±0.010 

0.0564 ± 0.0045 

17.208 ± 0.023 

8.37 ± 0.21 

0.585 ± 0.056 

0.957 ± 0.081 

JO106-1000 

16970 ± 260 

6.10 ±0.05 

0.191 ±0.010 

0.0642 ± 0.0051 

19.595 ± 0.023 

7.45 ± 0.20 

2.691 ± 0.253 

0.497 ±0.168 

JOl12+1835 

9740 ± 140 

5.77 ± 0.05 

0.160 ±0.010 

0.0863 ± 0.0080 

17.110 ±0.016 

8.01 ± 0.25 

0.662 ± 0.077 

1.822 ±0.196 

J0152+0749 

10800 ± 180 

5.93 ± 0.05 

0.168 ±0.010 

0.0736 ± 0.0061 

18.033 ± 0.009 

8.08 ± 0.22 

0.980 ±0.102 

1.407 ±0.137 

(This table is 

available in its entirety in a machine-readable form 

in the online journal. A portion is shown here for guidance regarding its form and content.) 

Table 5 

ELM WD Binary Parameters 



SDSS 

P 

(days) 

K 

(km s“^) 

Mass Function 
(M©) 

M2 

(Mo) 

^2,i=60° 

(Mo) 

'T’merge 

(Gyr) 

a 

(Re) 

log/l 

J0022+0031 

0.49135 ± 0.02540 

80.8 ± 1.3 

0.027 ± 0.003 

> 0.23 ± 0.02 

0.28 ± 0.02 


2.32 ±0.11 

-22.79 

J0022-1014 

0.07989 ± 0.00300 

145.6 ± 5.6 

0.026 ± 0.004 

5: 0.21 ± 0.02 

0.25 ± 0.02 

C 0.616 

0.65 ± 0.03 

-22.55 

J0056-0611 

0.04338 ± 0.00002 

376.9 ± 2.4 

0.241 ± 0.005 

> 0.46 ± 0.02 

0.61 ± 0.02 

C 0.120 

0.45 ± 0.01 

-22.06 

JO106-1000 

0.02715 ±0.00002 

395.2 ± 3.6 

0.174 ±0.005 

> 0.39 ± 0.02 

0.51 ±0.02 

C 0.036 

0.32 ± 0.01 

-22.61 

JOl12+1835 

0.14698 ±0.00003 

295.3 ± 2.0 

0.392 ± 0.008 

> 0.62 ± 0.02 

0.85 ± 0.03 

^ 2.682 

1.08 ± 0.01 

-22.40 

J0152+0749 

0.32288 ± 0.00014 

217.0 ±2.0 

0.342 ± 0.010 

^ 0.57 ± 0.02 

0.78 ± 0.03 


1.79 ± 0.03 

-22.80 


(This table is available in its entirety in a maehine-readable form in the online journal. A portion is shown here for guidanee regarding its form and eontent.) 


Table 6 

Radial Velocity Measurements 


Object 

HJD 

(days - 2455000) 

^helio 
(km s“') 

J03084-5140 

1634.76726 

19.2 ± 13.6 


1716.63827 

-145.6 ± 19.9 


1716.64861 

-132.3 ± 22.6 


1716.65891 

-154.4 ± 22.8 


1716.66922 

-134.4 ± 26.4 


1946.81093 

9.7 ± 12.1 


(This table is available in its entirety in a machine-readable form in the online 
journal. A portion is shown here for guidance regarding its form and content.) 
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